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Wave  Equation  Solution  for  Electromagnetic  Field  Incident  on 

Planar  Conducting  Arrays 

Peter  Loschialpo 

Abstract-A  new  wave  equation  solution  is  presented  for 
analyzing  planar  periodic  patch  or  aperture  arrays  having 
complex  conductivity  with  an  arbitrary  distribution.  The 
approach  extends  prior  developments  in  photonic  band-gap 
theory  by  including  the  conductivity  mechanism.  The  wave 
equation  is  written  in  the  form  Schrodinger ' s  equation. 
Solutions  are  obtained  in  the  spectral  domain  for  the 
radiated  field  in  terms  of  a  s^lmmation  of  plane  waves.  The 
method  provides  valuable  physical  insights  not  readily 
obtainable  using  conventional  solutions  for  periodic  arrays. 
Incident  and  radiated  fields  are  related  in  a  single 
mathematical  expression.  This  expression  shows  how  each 
radiated  field  mode  is  coupled  to  the  others  via  mutual 
capacitive  and  inductive  terms.  The  theory  may  be  readily 
extended  to  multiple  layers. 


I.  INTRODUCTION 

Due  to  their  wealth  of  applications,  two-dimensional, 
planar,  periodic  conducting  arrays  have  received  considerable 
attention  [1-9] .  They  have  recently  been  used  in  the  design 
Manuscript  approved  May  18, 1999. 
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of  reflectors  for  lightweight  spacecraft  antennas [8].  The 
design  exploits  resonances  in  order  to  separate  multiple, 
closely  separated  frequency  bands .  Accurate  predictions  in 
the  presence  of  grating  lobes  is  necessary  for  design  of 
frequency  scanning  antennas  [9]  which  are  designed  to 
efficiently  couple  energy  into  the  first  order  grating  lobe. 
These  periodic  structures  may  also  be  used  as  filters  or 
polarizers.  They  can  be  applied  over  a  wide  range  of  the 
electromagnetic  spectmm,  from  the  visible  to  the  microwave 
band. 


The  typical  approach  to  analyzing  such  structures, 
reviewed  in  detail  by  Mittra  [2]  and  Chan  [3],,  uses  the 
electric  field  integral  equation  (EFIE)  relating  the  incident 
electric  field  to  the  current  density  in  the  conducting 
patches.  Using  method  of  moments  techniques,  current  density 
is  expanded  as  a  sum  of  basis  fimctions  and  it's  coefficients 
are  found  by  the  solution  of  a  matrix  equation.  The  field 
radiated  by  the  patches  is  subsequently  calculated  from  the 
current  distribution.  In  a  similar  approach,  Usoff  [6]  and 
Munk  [7]  treated  the  current  array  as  a  sum  of  infinitesimal 
current  filaments .  These  current  filaments  are  written  in  a 
matrix  equation  as  a  function  of  the  voltage  induced  on  the 
filaments  by  the  incident  electric  field.  After  a  matrix 
inversion  to  solve  for  the  currents,  the  electric  field 
radiated  by  the  array  is  calculated  from  the  current. 
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Conductivity  is  assiomed  to  be  constant  within  a  patch  for 
these  approaches. 

A  different  method  is  described  which  offers 
potentially  valuable  new  physical  insights  as  well  as 
accurate  predictions.  In  a  similar  mcinner  used  to  investigate 
photonic  devices  [10-11],  the  wave  equation  is  written  in  a 
form  analogous  to  Schrddinger ' s  equation.  In  this  case,  the 
planar  array  has  a  periodically  varying  distribution  of 
conductive  patterns.  The  cortplex  conductivity,  or  surface 
admittance,  is  expanded  as  a  sum  of  exponential  basis 
functions.  Patterns  may  have  an  arbitrary  shape  and 
conductivity  distribution,  continuous  or  discrete,  provided 
enough  modes  are  included.  This  enhances  design  flexibility. 
Expressing  the  radiated  fields  as  a  sum  of  plane  waves,  a 
solution  is  derived  in  the  spectral  domain  consistent  with 
boundary  conditions  at  the  array  and  the  surrounding 
dielectric  regions.  The  solution  relates  the  incident  and 
radiated  fields  in  a  single  mathematical  expression,  with  the 
field  modes  coupled  to  each  other  by  terms  which  may  be 
identified  as  capacitive  and  inductive.  This  enables  a  direct 
examination  of  the  strength  of  the  interaction  between  field 
modes.  It  can  enable  a  designer  to  more  efficiently  couple 
more  energy  into  a  particular  mode  of  interest.  The  theory  is 
extendible  to  multiple  layers  in  a  straightforward  fashion, 
using  the  generalized  scattering  matrix  approach. 
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The  method  of  analysis  is  discussed  in  detail  in  Section 
II  for  a  conducting  planar  array  imbedded  in  a  multilayered 
dielectric  medium.  A  plane  wave  is  incident  to  the  dielectric 
stack  at  an  arbitrary  direction.  The  derived  matrix  equation 
can  be  solved  by  matrix  inversion  for  the  radiated  field 
modes.  These  are  then  used  to  calculate  transmission  and 
reflection  due  to  the  dielectric  stack  containing  the  array. 
Validation  measurements,  made  for  array  patterns  on  a  single 
layer  dielectric  are  discussed  in  section  III. 


II.  THEORETICAL  ANALYSIS 


The  analysis  of  the  interaction  of  a  planar  periodic 
array  of  conducting  shapes  with  an  incident  electromagnetic 
wave  is  in  investigated  through  a  wave  equation  in  the  form 
of  Schrodinger ' s  equation,  used  in  the  development  of 
photonic  band-gap  theory.  The  complex  admittance  distribution 
leads  to  a  matrix  solution  which  is  more  complicated  than 
that  for  the  photonic  band-gap  theory.  Field  modes  are 
coupled  to  on  another  by  capacitive  and  inductive  terms 
containing  the  Fourier  transform  coefficients  of  the 
admittance. 

Begin  with  a  thin,  planar  array  of  conducting  elements. 
The  elements  lie  in  the  plane  and  repeat  periodically 

with  distances  and  .  The  repeating  elements  may  have  any 
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shape  and  admittance  distribution  (see  Fig.  la) .  Figure  lb 
shows  this  array  imbedded  in  a  dielectric  slab  ”m”  having  a 
uniform  scalar  permittivity,  e„.  Slab  "m"  is  one  of  a  series 
of  M-1  slabs.  Free  space  surroxinds  multilayered  dielectric  on 
both  sides.  Subregions  of  slab  "m"  on  the  left  and  right  side 
of  the  conducting  array  are  identified  as  "a"  and  ”b", 
respectively . 

The  wave  equation  for  the  electric  field  in  any  of  the 
dielectric  regions  can  be  written  as 


-V^E  +  V(V  •  E)  -  ioyfij  =  <0 . 


(1) 


In  the  above  equation,  (O  is  angular  frequency,  is  the 

free-space  magnetic  permeeibility,  J  is  the  current  density, 

and  E  is  the  total  vector  electric  field.  MKS  units  are  used. 
Harmonic  time  dependence  is  assxomed  for  all  fields  and 

currents,  so  E(r,t)  =  E{r)e~^^^ .  Because  the  conducting  elements 
are  infinitesimally  thin,  the  z  corrponent  of  current  density, 
Jj,  is  zero.  The  x  and  y  components  are  given  by  Ohm's  Law 

written  in  the  form: 


Rsix,y) 


(2) 


Here  S(z-Zl)  is  the  dirac  delta  function  and  conductivity  can 
be  identified  as  5{z- Z[^)/  Rs(x,y) .  R^(x,y)  is  the  surface 
impedance  in  the  z  =  Zi  plane  of  the  array.  Note  the  similarity 
of  Eq.s  (1-2)  with  the  wave  equation  used  in  photonic  band- 
gap  theory  [10-11] .  The  periodically  varying  permittivity 
distribution  is  analogous  to  the  conductivity  in  Eq. (2) . 

The  electric  field  is  next  examined  in  detail  for 
dielectric  region  ” m" .  Subscripts  " m "  which  apply  to  the 
field  components  in  this  region  are  not  written,  though  are 
implicitly  understood.  Each  vector  component  of  the  field  may 
be  written  as  the  sum  of  an  "incident"  field,  E,.(r),  and  a 

"radiated"  field,  due  to  the  array  currents. 

£,(r)  =  E,(r)  +  Wr\r)  (  /  =  1,  2,  3  )  (3) 

Subscripts  ”  a”  and  "  b "  apply  to  the  two  subregions  of 
dielectric  slab  "  m "  on  either  side  of  the  array.  The 
subscript  i  =1,  2,  3  refers  to  the  three  vector  components 
jc ,  y ,  and  z  . 

The  incident  field  is  assumed  to  be  a  plane  wave.  It  is 
written  as 

E/(r)  =  .  ( 4 ) 
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In  Eq.  (4) ,  kj  is  the  component  of  the  wave  propagation 
vector  parallel  to  the  array.  It  is  given  by  kj  =  kix  +  k2y .  The 
component  perpendicular  to  the  array  is  given  by 

=  -^2  >  'Where  fil  =  e^{(o/cf  and  the  free  space  speed 

of  light  is  given  by  c . 

The  form  for  the  radiated  field  components  for  the  two 
subregions  "  a  ”  and  ”  b  ”  is  foimd  by  noting  previous  work 
[7] .  Munk  showed  that  the  field  radiated  by  a  periodic  array 
of  linear  current  elements  can  be  written  as  a  sum  of  plane 
waves.  The  same  is  true  for  any  current  distribution  which 
Ccin  be  expressed  as  a  sum  of  linear  current  elements .  The  x 
and  y  coitponents  of  the  radiated  field  are  then  written  as; 


G(I,2)- 


+i<:gj(z+Zj,) 


(5a) 

(5b) 


G  is  the  reciprocal  lattice  vector  of  the  array.  It  is 
given  by  G^GiX-h-Giy  and  its  conponents  are  ^x,y' 

where  ni2  =  0/  summation  over  G  is  taken  over 

both  components  Gj  and  G2  .  The  component  of  the  propagation 
vector  perpendicular  to  the  array  for  these  radiated  modes 
is : 
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(6) 


The  incident  field  coefficients  in  Eq.  (4)  are  related 
by  the  reflection  coefficients  at  the  z=z„  boundary.  This 

allows  Eq.  (4)  to  be  written  in  terms  of  the  coefficients  for 
the  negative  traveling  wave  alone. 


^1,2  (^)  ~  (^1,2- ■^0(1,2)  ^®-^(2,l)-^ 


+/^o,3  (^+  Zz.  )  \  ikj^  •r 


(7) 


+i*o^(z+Zj,-2^) 


(8) 


^(1,2)-  reflection  coefficients  (copolarization 
and  cross  polarization)  at  the  z„  dielectric  boiondary  for  the 
X  and  y  incident  field  components.  The  coefficients,  6^2-  in 
Eq.  (7)  are  the  input  variables  to  the  final  solution. 

Equations  similar  to  (7)  and  (8)  apply  to  the  field  in 
siibregion  "b"  radiated  by  the  array. 


U/*  /p-j  _  V*  C  !  »  U/fr 


(9) 


S  (?)  =  €  +  R  2Zm) 

^G(I,2) ^ 


C(l,2)- 


(10) 
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R  and  R-  are  the  reflection  coefficients  at  the  z  =  z„ 

G(l,2)-  C®- 

boiKidary  which  relate  the  negative  to  the  positive  traveling 
inodes . 


Equations  analogous  to  the  radiated  field  in  subregion 
"fe"  also  apply  to  the  field  in  siabregion  "a"  radiated  by  the 
array.  In  this  case,  and  R^^^  are  the  reflection 

coefficients  at  the  t  =  z„+i  boundary  which  relate  the  positive 
to  the  negative  traveling  radiated  modes. 


^1.2  (^)  ■"  2^('^C(1.2)+^C(1.2)  ^G®+'^G(2,I)+^ 


+  ^G(1,2).^ 


(11) 

(12) 


The  above  reflection  coefficients  for  the  X  and  y  field 
conponents  of  the  incident  and  radiated  fields  are  related  to 
the  well-known  reflection  coefficients  for  a  multi-layer 
dielectric  parallel  and  perpendicular  to  the  plane  of 
incidence  by: 


R 


G(l,2)± 


KAK-G^y+K,(k,-G,)^ 

{k,-Gy+(k,-G,f 


(13a) 


(13b) 
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on 


Coefficients  for  the  radiated  fields  'Ff  and  Y- 

G(l,2)+  C(l,2)- 

either  side  of  the  array  are  related  by  using  the  boundary 
condition  that  the  transverse  electric  field  is  continuous 
across  the  lattice. 


Y*  =n-  T-  +y-  Y? 

G(1.2)-  C(l,2)  C(I,2)+  '^0(2,1)  C{2,1)+ 


(14) 


The  dielectric  internal  reflection  factors  11-  and  y.  are 

C(l,2)  ^G(2,l) 

for  notational  convenience  and  are  defined  below  in  Eq.s 
(15a-b) . 


^0(1,2) 


(15a) 


'^C(l,2) 


S-  (z,)R-  (7 'SR 


G(h2)  ^ 
2 
G®- 


Aia^. 


(15b) 


The  constant  is  the  phase  shift  for  a  given 

mode  away  from  the  array.  Equations  (9),  (11),  and  (14)  allow 

the  field  radiated  by  the  array  to  be  expressed  by  a  single 
pair  of  field  amplitude  coefficients  ,  and  ,  in 

region  "m"  on  either  side  of  the  array. 


G2+ 
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To  derive  a  matrix  equation  for  the  radiated  field 
amplitude  coefficients,  and  incident  and 

radiated  fields  given  by  Eq.s  (7),  (9),  (11),  and  (14)  are 

substituted  into  Eq.  (1) .  The  resulting  equation  is  then 
integrated  term-by-term  over  z  from  Z/,-A  to  z^  +  A,  where  A 
is  a  vanishingly  small  number.  The  nonvanishing  terms  due  to 
the  discontinuity  at  the  array  lead  to  an  equation  which  can 
be  solved  to  obtain  the  desired  radiated  field  amplitude 
coefficients.  The  method  is  outlined  below. 

The  integral  of  the  first  term  in  Eq.  (1)  is 


A 


^u(?) 


dz 


(16) 


After  some  algebra,  the  integral  is  written  in  the  form 


z^+A 


a  ,  xua 

C(1.2)+  G(2,l)  0(2,1)+^ 


(17) 


The  dielectric  reflection-based  coefficients  in  Eq.  (17)  are 
defined  as 


(18a) 


(18b) 
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The  second  term  of  the  wave  equation  is  evaluated  by- 
consideration  of  Gauss's  law,  V»E  =  p/e,  and  the  continuity 

.  T-i  dp 

equation,  V»7  =  — These  lead  to 

dt 


V*E  =  - 


i 

(oe 


^dx  dy  j 


(19) 


Current  density  components  are  related  to  the  transverse 
field  and  surface  resistivity  by  Eq.  (2) ,  Since  the  surface 
admittance  of  the  array,  l^(r,y)  =  l/jRj(r,y) ,  is  a  periodic 

fiinction,  it  can  be  escpanded  as  a  sum  of  exponential  basis 
functions . 


(20) 


Equation  (2)  is  substituted  into  Eq.  (19) .  Both  surface 
admittance  and  the  radiated  field  coirponents  are  expressed  as 
a  s\im  of  the  exponential  basis  functions.  The  result,  after 
integrating  across  the  array  over  z,  is 
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Y[V(V  .  E)]  &  =  —  XX  Ji-(  -  «u)K m  -  G')Ej_S„,(z,)) 

Zf^-A  1  2  G  C'  ^  ^ 

+((*,  -  G,')£j.  +  ft  - 

+ft,,-G,,,-G'j)[(ift-G,-G;)TJ^C^ft))  (21) 

+te  -  G,  -  G.o'Tf,, + (t,  -  G,  -  got;,^  ] 

g-i(G+G')*r 


In  a  similar  way,  the  integral  of  third  term  of  the  wave 
equation  is  foimd  by  substitution  of  Eq.  (2) , 


Zi+A 


Z,-A 


G  O' 


(22) 


Because  the  transverse  components  of  the  electric  field 
are  continuous  across  the  array,  the  integral  of  the  fourth 
term  of  the  wave  equation  is  zero. 

The  z -integrated  wave  equation  is  arrived  at  by  adding 
Eq.s  (17),  (21),  and  (22)  for  the  individual  terms. 

Sijmmations  in  the  resulting  equation  are  carried  out  over 
both  G  and  G' .  To  obtain  a  form  which  can  be  solved  for  the 
radiated  field  amplitudes,  the  resultant  equation  is  first 

multiplied  by  ,  then  integrated  over  the  x  -  y  plane; 

The  exponential  basis  functions  are  orthogonal,  allowing  for 
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one  of  the  summations  to  vanish.  A  general  equation  for  the 
radiated  electric  field  amplitudes  is  found. 


+((t,  -  G,)£i.  +  (t,  -  G,  )£,_  ]Y^ 

^(K2  -  C,.2)  J[(i(t,  -  G,)'Pa,.Ca,fe)) 


(23) 


mo 


C,3 


This  general  equation  can  be  put  in  a  matrix  form  and  solved. 
Considerable  simplification  is  obtained,  however,  in  cases 
where  the  coupling  between  the  x  and  y  components  of  the 
electric  field  can  be  neglected.  Then  the  number  of  matrix 
elements  is  reduced  by  a  factor  of  four,  enabling  much  faster 
solutions . 


For  the  remaining  analysis,  neglect  the  coupling  between 
the  X  and  y  field  components.  Equations  (15a-b)  and  Eq.s 
(18a-b)  reduce  to 

C(l,2)  o  (  \ 

(24) 

Y-  =0 
^G(l,2) 
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and 


^(1,2)  ^0(1.2)+^ 
Tf  =0 

C(1.2) 


-2/a, 


'"')+ng„,„0-« 


2/a, 


0(1,2)- 


■dm^ 


(25) 


Equation  (23)  simplifies  to 


(26) 


Ecjuation  (26)  represents  a  set  of  2N  equations,  where  N  is 
the  product  of  the  number  of  x  and  y  reciprocal  lattice 

vectors .  Their  solution  yields  the  modes  and  which 

comprise  the  x  and  y  coitponents  of  the  radiated  electric 
field.  An  examination  of  Eq.  (26)  shows  how  each  mode  is 
coupled  to  the  other  modes  through  a  capacitive  term, 

proportional  to  — — ,  and  an  inductive  term,  proportional  to 

CDEk- 

G.3 


The  equation  shows  how  the  aitplitude  of  a  particular 


mode  is  affected  by  parameters  such  as  dielectric 
reflections,  incident  wave  direction,  and  pattern  shape.  For 
nxomerical  solution,  Eq.  (26)  is  written  in  the  form  of  a 
matrix  eq[uation.  ^ 
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(27) 


A  +  foB  +  C)  •  X  =  ft)"D  +  F . 


For  either  the  x  or  y  field  components,  elements  of  the 
square  matrices  A,  B,  and  C  are  identified  as 


A... 

=  -^C- 

G.G' 

ifc- 

G,3 

B-  . 

=  TrS-  - 

G,G 

1 

G,G' 

sk 

^C.3 

*  C,i  G-G 


(28) 


Elements  of  the  column  matrices  D,  F,  and  X  are 


^C,3 

^G.3 


X- =4^;. 

G  G,»+ 


(29) 


Equation  (29)  can  be  solved  using  standard  matrix 
inversion  techniques.  The  column  matrix  X  contains  the 
amplitude  coefficients  for  the  electric  field  Floquet  modes 
radiated  by  the  conducting  array.  These  modes  are  plane 
waves.  They  are  linked  together  by  the  coefficients  for  the 
array  conductivity,  .  Together  with  the  incident  field 

added  in  phase,  and  including  the  reflections  at  dielectric 
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boundaries,  reflection  and  transmission  coefficients  are 
calculated  in  a  straightforward  manner. 


III.  EXPERIMENTAL  VALIDATION 

Accuracy  of  the  theoretical  analysis  has  been  verified 
by  a  series  of  experimental  measurements.  Transmission 
measurements  for  a  number  of  test  devices  exhibit  a  variety 
of  striking  characteristics,  including  resonances,  which  are 
well  reproduced  by  calculations. 

Results  of  ejqjerimental  validation  studies  are  presented 
for  four  test  panels.  For  three  panels  the  patterns  were 
highly  conductive.  Correspondingly,  the  surface  resistivity 
entered  in  the  computations  was  chosen  small  enough  to  have 
no  significant  effect  on  the  output.  The  test  panels  are 
approximately  square  with  a  length  of  30  cm  on  each  side. 
Measurements  were  performed  in  a  transmission  tunnel  which 
contains  two  antennas  for  transmission  and  reception  of 
microwave  radiation  inside  a  box  lined  with  absorber 
material.  A  baffle  following  the  transmitting  antenna 
collimates  the  incident  wave.  Figure  2  displays  a  schematic 
of  the  experimental  arrangement.  Incident  energy  is  normal  to 
the  array  pcinel  with  the  electric  field  polarized  in  the  x 
direction,  parallel  with  one  edge  of  the  square  patte3ms. 
Measurements  were  made  using  a  Hewlett  Packard  8 5 IOC  network 
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analyzer.  Effects  of  spurious  reflections,  such  as  those  off 
the  wall,  were  minimized  by  suitably  gating  the  signal  in  the 
time  domain. 

The  first  panel  was  fabricated  from  a  commercial 
photosensitive  printed  circuit  board.  After  etching,  a 
pattern  of  thin  copper  squares  remained  on  the  0.145  cm  thick 
dielectric  backing.  The  periodicity  of  the  pattern  is  1.0  cm 
in  each  direction  and  the  copper  square  size  is  0.9  cm. 
Electric  permittivity  (normalized  to  free  space)  of  the 
dielectric  was  measured  and  found  to  have  a  real  part  of  4.4 
and  an  imaginary  part  less  than  0.1.  Transmission  was 
measured  with  the  copper  square  pattern  facing  the  incident 
wave.  Matrix  calculations  included  129  modes  in  the  direction 
of  the  incident  field  and  17  modes  perpendicular  to  it  (2193 
unknowns) .  Resulting  curves  for  both  experiment  and 
conputation  are  shown  in  Fig.  3.  As  expected,  the  test  panel 
acts  as  a  low  pass  filter,  with  transmission  decreasing  with 
increasing  frequency.  Calculated  and  measured  transmission 
differ  by  at  most  0.05  over  the  entire  frequency  range.  The 
frequency  at  which  the  wavelength  of  a  plane  wave  propagating 
in  the  dielectric,  is  equal  to  the  array  period,  is 

noted  on  the  graph.  This  is  the  frequency  above  which  the 
lowest  order  secondary  modes,  =  (±1,0)  and  (0,  ±1)  ,  can 

propagate  in  the  dielectric,  as  can  be  seen  from  Eq.  (6) . 
Physically,  the  array  constrains  these  modes  to  have  a 
periodicity  which  is  a  multiple  of  in  the  ^  direction. 
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irt5)lying  the  x  coirponent  of  the  propagation  vector  is  2nld^. 

The  vector  sxom  of  the  x  component  and  the  z  coirponent  (must 
be  real)  of  the  propagation  vector  must  have  a  magnitude  of 
IniXjy.  This  only  is  possible  ifAo<</,.  In  the  literature, 

these  secondary  modes  are  known  as  a  grating  lobes .  Because 
of  the  total  internal  reflection  by  the  dielectric  slab,  the 
grating  lobes  are  trapped  inside  the  dielectric  until  the 
frequency  is  above  the  point  at  which  the  free  space 
wavelength,  is  less  than  d^.  The  (±1,0)  and  (0,±1) 

grating  lobes  then  propagate  away  from  the  panel  in  a 
nonspecular  direction.  In  general,  higher  order  modes 
propagate  when  as  given  by  Eq.  (6),  is  real  in  the  medium 

of  interest. 

A  second  test  sairple  consists  of  an  array  of  aluminum 
foil  squares  with  a  periodicity  of  2.0  cm  in  each  direction. 
The  array  is  taped  to  a  0.01  cm  thick  paper  backing.  Total 
dielectric  thickness,  tape  plus  paper,  is  approximately  0.02 
cm.  Lengths  of  the  squares'  edges  are  1.5  cm.  For 
conputations,  the  small  effect  of  the  paper  backing  was 
included  by  using  a  real  permittivity  of  5.1,  found  by 
earlier  measurements.  The  tape  was  neglected.  Calculations 
included  65  modes  in  the  direction  of  the  incident  field  and 
17  modes  perpendicular  to  it.  Calculated  and  measured  values 
of  transmission  are  graphed  in  Fig.  4  (solid  lines  and 
circles) .  The  panel  is  a  low  pass  filter  below  the  resonant 
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frequency,  CO  =  9 . 0x10^°  rad/s .  The  primary  mode,  (0,0),  is  the 
only  propagating  mode  up  to  this  frequency.  At  resonance,  the 
amplitude  of  the  primary  radiated  mode  peaks  with  an 
amplitude  and  phase  which  nearly  cancel  the  incident  wave.  An 
examination  of  the  field  modes  has  shown  that  is  due  to  the 
coupling  of  the  incident  field  to  the  evanescent  fields  which 
has  a  90°  phase  shift  at  this  frequency.  A  second  90°  phase 
shift  occurs  in  the  coupling  from  the  evanescent  fields  to 
the  radiated  field,  resulting  in  the  net  180°  phase  difference 
between  the  incident  and  radiated  primary  field.  A 
substantial  loss  in  transmission  results,  which  is  clear  in 
the  figure.  For  a  given  periodicity,  we  obsesrved  that  the 
resonant  frequency  is  dependent  on  pattern  dimensions.  The 
resonant  behavior  can  be  modeled  as  a  circuit  analog.  Figure 
4  (dashed  lines  and  x's)  also  shows  the  results  of 
measurements  and  calculations  for  the  same  geometry,  but  with 
the  conductive  patterns  cut  from  a  sheet  having  a  surface 
resistivity  of  53  Qp.  As  expected  from  circuit 

considerations,  the  resistivity  dampens  the  resonant  peak, 
but  otherwise  the  shape  of  the  curves  are  similar.  Above  ft)  = 
9.4x10^°  rad/s,  and  the  modes  (±1,0)  and  (0,±1)  begin 

to  propagate  away  from  the  panel.  This  follows  from  the 
arguments  given  in  the  previous  paragraph.  These  modes 
propagate  in  a  nonspecular  direction  with  a  polar  angle  given 
by  cosd  =  k-^l and  therefore  to  not  contribute  directly  to  the 

transmission  coefficient.  The  excellent  agreement  with 
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measured  data  is  also  evident  in  Fig.  4.  Resonances  in  the 
transmission  spectrum  have  historically  been  referred  to  as 
Wood's  anomalies,  which  were  associated  with  the  onset  of 
propagation  of  grating  lobes  in  slot  arrays  for  the  free 
space  case  [12].  The  resonances  here  are  analogous. 

The  fourth  panel  consists  of  an  a  array  on  1.5  cm 
aluminiim  squares  with  a  2.0  cm  periodicity  like  in  the 
previous  figure.  In  this  case,  it  is  backed  by  a  0.54  cm 
thick  Plexiglas  sheet.  Calculations  included  65  modes  in  the 
direction  of  the  incident  field  and  17  modes  perpendicular  to 
it.  Calculated  and  measured  results  are  compared  in  Fig.  5. 
The  addition  of  a  dielectric  backing  has  resulted  in  two 
resonance  peaks,  one  at  7.0x10^°  and  the  other  at  8.8x10 
rad/s.  The  (±1,0)  and  (0,±1)  modes  are  trapped  by  the 
dielectric  in  the  frequency  range  from  6.2x10^°  {X^  =  d^)  to 

9.4x10^°  rad/s  =  Above  8.8x10^°  rad/s,  =  and  the 

(±1,±1)  modes  can  propagate  in  the  dielectric,  but  not  in 
free  space.  These  frequency  intervals  are  indicated  in  the 
figure.  Agreement  is  quite  good  between  the  measured  and 
calculated  results.  Observed  discrepancies  have  a  number  of 
possible  causes.  Calculations  assime  the  planar  array  is 
infinite  in  extent.  Yet,  the  baffle  which  limits  the  beam 
size  is  20  cm  in  diameter,  equivalent  to  10  array  periods  for 
this  panel.  Departures  of  the  incident  beam  wave  front  from 
the  assimied  plane  wave  can  also  be  significant. 
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IV.  CONCLUSION 


A  new  method  for  analyzing  the  electric  field  scattered 
by  a  planar  periodic  array  of  conductive  elements  in  a 
stratified  dielectric  has  been  developed.  The  analogy  of 
these  structures  with  photonic  devices  was  exploited  in 
formulating  the  technique.  Coupling  between  radiated  field 
modes  as  well  as  with  the  incident  field  is  given  by  a  single 
equation  which  can  be  solved  by  matrix  inversion  techniques . 
Since  mode  coupling  determines  overall  performance,  this 
approach  may  lead  to  better  understanding  of  scattering  by 
these  structures.  Periodically  varying  conducting  elements 
are  expressed  as  a  sum  of  exponential  basis  fiinctions, 
allowing  for  arbitrarily  shaped  complex  impedance 
distributions  with  a  continuous  or  discontinuous  profile  to 
be  used.  This  enables  increased  design  flexibility. 

Theoretical  predictions  were  made  for  a  variety  of 
configurations,  from  a  conputer  program  based  upon  this 
method.  Accuracy  of  the  method  was  demonstrated  by  comparing 
experimental  transmission  measurements  with  the  model's 
predictions.  Results  were  in  excellent  agreement  with 
calculations .  Frequencies  at  which  grating  lobes  appear  were 
noted.  This  method  can  be  extended  to  multilayered  periodic 
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structures,  using  the  generalized  scattering  matrix  approach, 
in  a  straightforward  manner. 
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FIGURE  CAPTIONS 


1.  (a)  Thin,  conducting  array  in  the  x  -  y  plane  with 
periodicity  and  .  (b)  The  planar  array  at  z  =  Zl 
imbedded  in  a  multilayered  dielectric  medium  ( z,  to  z^,). 

2.  Experimental  configuration  used  to  measure  the  transmission 
coefficients . 

3.  Conparison  of  calculated  with  measured  values  of 
transmission  through  an  array  of  copper  squares  having  a 
1,0  cm  periodicity  backed  by  a  0,145  cm  thick  dielectric 
sheet , 

4.  Comparison  of  measured  and  calculated  transmission  for  two 

different  square  arrays  (aluminiam  and  and  53  each 

with  a  2,0  cm  periodicity.  Calculated  and  measured  values 
for  the  periodic  Al  array  are  represented  by  solid  lines 
and  open  circles.  The  resonant  dip  at  90x10®  rad/s  is  due 
to  the  coupling  between  the  incident  and  evanescent  fields 
which  are  90°  out  of  phase  at  this  point.  For  the  periodic 
resistive  array,  calculated  and  measured  transmission  are 
represented  by  dashed  lines  and  x's, 

5.  Cortparison  of  calculated  and  measured  transmission  through 
an  array  of  alumin\am  squares  having  a  2,0  cm  periodicity. 
The  array  is  backed  by  a  Plexiglas  sheet  having  a  thickness 


of  0,54  cm. 
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Figure  4 
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Figure  5 
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